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An effective description for spherical nanoparticles in a fluid of point particles is 
presented. The points inside the nanoparticles and the point particles are assumed 
to interact via spherically symmetric additive pair potentials, while the distribution 
of points inside the nanoparticles is taken to be spherically symmetric and smooth. 
The resulting effective pair interactions between a nanoparticle and a point parti- 
cle, as well as between two nanoparticles, are then given by spherically symmetric 
potentials. If overlap between particles is allowed, the effective potential generally 
has non- analytic points, but for each effective potential the expressions for different 
overlapping cases can be written in terms of one analytic auxiliary potential. Ef- 
fective potentials for hollow nanoparticles (appropriate e.g. for buckyballs) are also 
considered, and shown to be related to those for solid nanoparticles. Finally, explicit 
expressions are given for the effective potentials derived from basic pair potentials of 
power law and exponential form, as well as from the commonly used London- Van der 
Waals, Morse, Buckingham, and Lennard-Jones potential. The applicability of the 
latter is demonstrated by comparison with an atomic description of nanoparticles 
with an internal face centered cubic structure. 

PACS numbers: 62.23.Eg, 36.40.-c, 02.30.Mv 



I. INTRODUCTION 

Nanoparticles,^^ quantum dots, 4 colloidal suspensions,^ and globular proteins^ are ex- 
amples of physical systems in which small nanometer or micron-sized clusters of particles are 
suspended in a fluid. Such systems have applications ranging from material coatings to drug 
deliveryPSl For colloidal systems, collective behavior has been the focus of much research)*^ 
while nanoclusters are often studied as isolated objects " 1 * 12 * 13 * 14 * 1 -^ despite interesting collec- 
tive phenomena such as the increased heat conductance in dilute nanoparticle suspensions^ 
and self-assemblyP 

To study the collective properties of nanoparticles in suspension, one would expect that 
a detailed description of the internal structure of the clusters is not necessary, especially if 
the nanoparticles are more or less solid. On the other hand, a description in terms of hard 
spheres would probably be too crude for nanoparticles since typical atomic interaction ranges 
are on the order of Angstroms. The aim of this paper is to give a general effective description 
of nanoparticles which retains a level of detail beyond the hard sphere model and which is 
intended to be used in the study of the collective behavior of nanoparticles, either numerically 
or analytically. The starting point of the description is to assume that each nanoparticle 
is composed of particles with fixed relative positions, interacting with the point particles in 
the fluid and their counterparts in other nanoparticles through spherically symmetric pair 
potentials. It is furthermore assumed that the nanoparticles may be modeled as spheres 
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with a smooth spherically symmetric density of constituents, which can be viewed as a 
smoothing procedure for the interactions. In particular, solid and hollow spheres of uniform 
density are considered in detail, since these are suitable for describing solid nanoclusters and 
buckyballs (or similar structures), respectively. The spherical smoothing procedure results 
in spherically symmetric effective interaction potentials for nanoparticles and point particles, 
and consequently leads to a description of a nanoparticle as a single particle instead of as a 
collection of particles. 

Similar approaches to the problem of constructing effective potentials have been used 
before, but only for specific cases JUMHHSHSni] The current paper is devoted to the general 
method of deriving effective pair potentials for nanoparticles from the basic pair potential 
of their constituents. The possibility of overlapping and embedded particles is specifically 
treated as well. 

The paper is structured as follows. In Sec. [TT| the general smoothing procedure is ex- 
plained. Properties of the resulting effective potentials are explored in Sec. Ill with special 
consideration for the difference between non-overlapping and overlapping particles, which re- 
sults in a reformulation of the non-analytic effective potentials in terms of analytic auxiliary 
potentials. In Sec. IV, the formalism is extended to include hollow nanoparticles. For uni- 
form solid and hollow nanoparticle structures, explicit effective potentials for a nanoparticle 
and a point particle and for different nanoparticles are worked out in Sec. [Vjfor the London- 
van der Waals potential, the exponential potential, the Morse potential, the (modified) 
Buckingham potential, and the Lennard- Jones potential. Section VI addresses the applica- 
bility of the effective potentials by comparison with an atom-based nanoparticle model. A 
discussion in Sec. |VII| concludes the paper. 



II. SMOOTHING PROCEDURE FOR NANOPARTICLE POTENTIALS 

Consider a classical system of point particles, representing a fluid, and spherical clusters 
called nanoparticles. While in reality, a nanoparticle is a cluster of a number of atoms, here 
each nanoparticle will be modeled by a smooth internal density profile p(x) that depends 
on the distance x from the center of the nanoparticle only and which is strictly zero for 
x > s, where s is the radius of the spherical nanoparticle. This approximation is motivated 
by the idea that for spherical nanoparticles, the inhomogeneities due to the discreteness 
of the atoms inside the nanoparticles should only have a small influence on the effective 
nanoparticle potentials. Given a density profile p(x), one can make contact with the picture 
of a nanoparticle as a cluster of distinct atoms by interpreting M = /g dx p(x) as the total 
number of atoms inside the nanoparticle, where x — |x|, and B s denotes that the integration 
over x is over the volume of a ball of radius s around zero. 

To further illustrate that it is reasonable to smooth out the internal density, consider the 
idealized case that the atoms composing the nanoparticle are arranged in a face-centered- 
cubic (fee) lattice — the crystal structure of e.g. aluminium, silver, gold, and platinumP^ — 
with one of the atoms in the center. The true density inside the nanoparticle is then a sum 
of delta functions, but this can be coarse-grained by taking a spherical shell of radius x with 
a width 8x, counting the number of atoms in the shell, and dividing by the volume of the 
shell. The result of such coarse-graining is shown in Fig. [T] for a lattice with mean number 
density p = 1 and for two values of the coarse-graining width, 5x = 3/4 and Sx = 3/2. 
The coarse-grained density around a single atom in an fee crystal is seen to be reasonably 
constant except near the central atom (with the positive and negative deviations from the 
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FIG. 1: Coarse-grained radial density profile of the fee lattice of mean density p = 1 as a function 
of the distance from a central atom. The circles correspond to a coarse-graining width of 5x = 3/4, 
the squares corresponds to dx = 3/2 (the points are connected to guide the eye). The horizontal 
line indicates the mean number density. 



mean density averaging out for larger Sx), so that to first order the density may be replaced 



by a constant. This highly idealized nanoparticle structure will be used again in Sec. |VI| to 
get an idea of the accuracy of the effective potentials. 

Let 4> pn {f) denote the basic pair potential between a point of a nanoparticle and a point 
particle in the fluid, where r is the distance between them. This potential will be assumed 
to be analytic for r > but may diverge as r — > 0. The effective point-nanoparticle pair 
potential V pn is then given by 



V pn (r)= / dxp(x)0 pn (|r-x|), (1) 

•IBs 

where the subscript pn denotes that this is a point particle-nanoparticle potential and r is 
the distance vector between the point particle and the center of the nanoparticle. Because 
of the spherical symmetry of the density profile and the pair potential pn , the effective 
potential does not depend on the direction of r, only on its magnitude r = |r|. 

Analogously, the effective inter-nanoparticle potential V nn for two nanoparticles with in- 
ternal density profiles p\ and p2, radii si and S2, and whose points interact through a pair 
potential nn , is given by 

V m (r) = / dx / dy px(x)p 2 (y) </> nn (|r - x - y|), (2) 
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The potential nn will also be assumed to be analytic for r > 0. Throughout this paper, pn 
and nn will be referred to as the basic pair potentials, while V pn and V nn are the effective 
potentials. 

To arrive at more concrete expressions for the effective potentials, it will be assumed that 
the internal density profile of the nanoparticles is analytic, so that it may be written as a 
Taylor series, 

oo 

p(x) = 6(s - x) E a i x \ ( 3 ) 

i=0 
i even 

where B is the Heaviside step function. In Eq. (|3]), odd powers of x were omitted since 
they lead to non- analytic behavior at x — 0. The potentials for a nanoparticle and a point 
particle, and for two nanoparticles, respectively, that would result from internal densities of 
monomial form 0(s — x)x % are denoted by 

Vi(r) = I dxx 4 pn (|r + x|), (4) 
Jb s 

Vij{r) = I dx/ dya;V> nn (|r + x-y|). (5) 

Here, and below, the dependence of Vi and Vij on s and si and $2 will not be denoted 
explicitly. In terms of the potentials Vi and V^, the effective point-nanoparticle and inter- 
nanoparticle potentials are given by 

oo 

^pn(r) = £ a^(r) (6) 

i=0 
i even 
oo oo 

Kn(r) = E E «W(r) (7) 

j=0 j=0 
i even j even 

where = 0(si— x) X)j ^j^ 1 an d ^(^c) = 0(^2— x) X),- ^i^' 7 are the internal density profiles 
of two interacting nanoparticles. While often only the first term i = j = will suffice, the 
formalism will be developed for general i and j, since this is not any more difficult. 

The three-dimensional and six-dimensional integrals in Eqs. Q and ^ for the effective 
potentials make further manipulations cumbersome. However, due to the spherically sym- 
metry of the basic pair potentials, these multi-dimensional integrals can be rewritten as 
integrals over a single variable. 

To convert Eq. (El) to a single integral, one goes over to spherical coordinates x = 
(xsin^cos^xsin^sin^a; cos#), integrates over tp and then performs a change of inte- 
gration variable from 9 to y = [x 2 sin 2 9 + (r — x cos 9) 2 } 1 ^ 2 , which yields 

v i{r) = — dx dyx l+1 y<b pn (y), 



r JO J\r-x\ 

Reversing the order of the x and y integrals and using that i is even leads to 



VAr) 277 



(i + 2)r 



r+S dy [s^ 2 -(r-yY +2 } y <P VR { y ) 

\r-s\ 



+Q(s - r) / dy[(r + yY +2 - (r - y) 4+2 ]# pn (y) 



(8) 



5 



Defining a kernel 



KAx,s) 



2tt 
i + 2 



J+2 



x i+2 )e(s- \x\), 



one can write the right hand side of Eq. ([8| in the concise form 



VAr) 



dyKi(r-y, s) y <p pa (\y\) 



(9) 



(10) 



at least for r > s. That Eq. (10) also holds for r < s (with the same expression for Ki) is 
seen by writing the second term in Eq. ([8| as 



dy[W +2 - (r - y y+ 2 } - {s^ 2 - (r + yy +2 }]y<P P M 
^ dy[s i+2 -{r-y) i+2 )y4> m {\y\)- 



s+r 



Combining this with the first term in Eq. (|8j) leads again to Eq. (10). Note that for the 
special case of i = 0, to be used below, the kernel takes the form 



K (x, s) = 7r(s 2 — x 2 ) 0(s — \x\). 



(11) 



For the effective inter-nanoparticle potential V^-, one can use that the potential energy 
of two nanoparticles is equivalent to the potential energy of a particle and a nanoparticle of 
which the points interact via a point-nanoparticle potential Vj, i.e., 



dy K i (r-y,s 1 )yV j (\y\), 



where in Vj, one should replace s by s 2 , and pn by nn . Combining this with Eq. (10), and 
using that Kj(x, s 2 ) is even in x, one obtains 



Vij{r) = - JdydzK i (r-y,s 1 )K j (y-z,s 2 )z(j) I1Il (\z\), 



or 



dyK i:j (r - y, s u s 2 ) y (f> n n(\y\) , 



with the kernel Ku given by 



Kij(x, Si, s 2 ) = Jdy Ki(x - y, si) Kj(y, s 2 ). 



(12) 



(13) 



(14) 



The integral in this expression is further evaluated in the Appendix, where it is shown 
that Kij is a piecewise polynomial function of degree % + j '• + 5 which has a finite support 
M < Si + s 2 , and non-analytic points at x = ±|si — s 2 \. For the special case i = j = 
which will be used below, one finds from Eqs. (A3) and (A4), and after some rewriting, 

fjGD - \d\) 3 {d 2 + 3D\d\ + D 2 - 5x 2 ) if |x| < \d\ 

l\ {){l (r..^.s. 2 ) { g( jD _| x |)3( x 2 + 3jD | a; | +jD 2_ 5rf 2) if |d| < | x | < £) (15) 

if b| > D, 
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where 



D = si + s 2 

d = si — s 2 - (16) 



Because the kernels Ki and Kij are piecewise polynomials, the integrals in Eqs. (10) 



and (13) can be performed analytically for many functional forms of pn and nn , such as 
power law and exponential forms (see Sec. [V]), which are the basis of many commonly used 
empirical pair potentials. 

III. AUXILIARY POTENTIALS 



Although not evident from Eqs. (10) and (13), the non-analytic points of the kernels and 



of the basic pair potential cause the effective potentials to have different functional forms 
depending on whether there is overlap. Different overlapping cases can occur: A point 
particle and a nanoparticle can either overlap (for r < s) or not overlap (for r > s), while 
two nanoparticles can have no overlap, which requires r > s\ + s 2 = D, or partially overlap, 
or the smallest nanoparticle can be completely embedded in the larger, which occurs when 
r < \s\ — s 2 \ = \d\. The different forms of the effective potentials for these different cases 
can be linked by introducing auxiliary potentials. 

The following symmetrization operations on functions / are useful in denoting the rela- 
tions between effective and auxiliary potentials 

f([x]) = f(x) — f(—x) "antisymmetrization" 
f((x)) = f(x) + f(—x) "symmetrization." 

These operations are also useful for functions with multiple arguments, e.g., 

/(M>2/) = ffav) - f(-x,y) 

f( x ,(y)) = f(x,y) + f(x,-y) 

/(NjM) = f(x,y) - f(-x,y) - f(x,-y) + f(-x,-y) 
f(i x ,y]) = f( x ,y) - f(-x,-y). 

Note that in the last example, a single antisymmetrization was performed which involved 
both arguments. 

The expressions of the effective potentials Vi and V^- in terms of auxiliary potentials 
(whose derivations will follow) are given by 

{Aij((r), [si],s 2 ) if r < \d\ and si < s 2 

A ij{{r), si, [s 2 ]) if r < d and si > s 2 , . 

^•((r), 8l , s 2 ) - Ai^r, ( Sl , -s 2 )) if |d| < r < D ^ 
A ij( r , [si], N) if r > D, 

in which the auxiliary potentials are defined as 

1 rr+s 

Ai(r,s) = - dyKi(r-y,s)y4> pn (y) (19) 
1 J 

\ rr+si+S2 

Aij(r, s h s 2 ) = - / dy Kij(r -y,s 1 ,s 2 )y(p nn (y), (20) 

/ J 
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where furthermore 



KAx, s) 



2tt 

2 + 2 

fX + Sl 



(s l+2 - X 



i+2\ 



/x+si 
dy Ki(x - y, si) i^(y, s 2 ). 
-82 



(21) 
(22) 



Note that is the analytic continuation of K iy while the quantity Kij(x, si, 8%) has the 
same functional form as the kernel forx<0,G?<|x|<.D (as it coincides with case 4 
in the appendix). In particular, for i = j = 0, one has from Eq. (15) 



K Q0 (x,si,s 2 ) 



7T 

30 



(si + s 2 + x) 3 {x 2 - 3six - 3s 2 x - As\ - 4s^ + 12sis 2 ). (23) 



The derivation of Eq. (17) goes as follows. Consider first the non-overlapping case r > s. 
In that case, the absolute value sign in the argument of pn may be dropped in Eq. ( 10 ), since 



r > s and r — y < s [cf. Eq. rt9j)] imply that y > 0. Thus, the effective point-nanoparticle 
potential can be written as 



Vi(r) 



- dy Ki(r — y,s)y pn 

(y) 

\ rr+s 

- / dy Ki(r -y,s)y<j) pD (y) 

T Jr—s 

1 rr+s _ 1 rO 

- / dy Ki(r -y,s)y (p pn (y) + - dy K^r -y,s)y 4> pa (y) 
r Jo r Jr-s 



Ai(r, s) - Ai(r, -s) 



(24) 



For the case r < s, the argument in the pn function in Eq. (10) needs to be —y for y < 0, 
giving 



Vi(r) 



X rr+s | rO 

/ dy Ki(r -y,s)y<f) pn (y) + - dy K^r - y, s) y <f) pn {-y) 

J f J T ~ S 

dyKi(-r-y, s)y<p pn (y), 



r Jo 
1 



r Jo 



r+s _ | 

dy i^j(r - y, s)y<f) pn (y) - - 

r Jo 



r—s 
s—r 



(25) 



where a change of integration variable from y to — y was carried out in the second integral, 
and it was used that Ki(y, s) is even in y. The first term on the right hand side of Eq. (25) 
is equal to Ai(r, s) in Eq. (19), while the second term equals r, s), so that 

Vi(r) = A(r, s) + M-r, s) = M(r), s). (26) 

Thus, although the effective potentials between a point particle and a nanoparticle have 
different forms for non-overlapping and overlapping situations [Eqs. (24) and (26), respec- 
tively], both can be written in terms of the auxiliary potential A4, and one obtains Eq. (17). 

A technical difficulty must be mentioned here, namely, that the integral defining the 
auxiliary potential in Eq. (19) may not converge, even when the linear combinations in 
Eq. (17) do. In such cases, one should strictly write the auxiliary potential as a sum of a 
regular and a diverging part by replacing the lower limit of the integral in Eq. ( 19 ) by 5 > 0, 
and expanding the result in 5. In the absence of overlap, Eq. (17) must yield a finite result, 
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FIG. 2: Subdivision of the integration domain in the derivation of the expression of the inter- 
nanoparticle effective potential Vy in terms of the auxiliary potential Aij. Assuming si > S2, three 
cases have been distinguished: (a) r > s\ + s 2 , (b) si — S2 < r < si + S2, and (c) r < si — s 2 . 



i.e., the diverging parts (negative powers of S and possibly logarithmic terms) must cancel, 
hence in that case it suffices to work with the regular part of the auxiliary potential. On the 



other hand, in case of overlap, it is possible that the divergent parts do not cancel in Eq. (26), 
resulting in infinite effective potentials. An independent criterion for whether an effective 
potential is infinite in overlapping cases can be constructed as follows. For a single particle 
inside a nanoparticle, the effective potential becomes infinite only if the divergence of the 
basic pair potential pn at the origin is too strong. In particular, if 0(r) oc r~ k for small r 
then the point-nanoparticle potential is infinite for k > 3, as is seen by considering a small 
sphere around the particle, giving an integral of the form f r<s dr 0(r) oc J^dr r 2 r~ k ~ 
which diverges for k > 3 in the limit 5 — > 0. This result extends to inter-nanoparticle 
potentials, which are also infinite if there is overlap and the potential nn diverges no slower 
than r~ 3 , i.e., the Vij(r) are finite for r < D provided 4> nn (r) diverges for small r slower 
than r~ 3 . Given this criterion, the divergent part of an auxiliary potential is not needed 
to determine whether the corresponding effective potential is infinite. Since the divergent 
parts are needed neither in overlapping nor in non-overlapping cases, below, only the regular 
parts of auxiliary potentials will be given. 



To derive Eq. (18) for the effective potentials between two nanoparticles, one starts by 



rewriting Eq. (12) to 

v ij( r ) = ~ d V dxK i (y,s 1 )K j (x,s 2 )(r-x-y)(f) mi (\r-x-y\). (27) 

T J— si J — S2 

In this formulation, the integration domain is a rectangle in the (x, y) plane and the integrand 
has a diagonal non-analytic line at x + y = r. This line may or may not cross the domain, 
which is what gives rise to non-analyticity and the difference between overlapping and non- 
overlapping effective potentials. 

Subdividing the domain into triangular regions without non-analyticities will result in 
expressions in terms of analytic subexpressions. The appropriate subdivisions of the inte- 
gration domain are shown in Fig. |2j where it was assumed that the radius si is larger than 
the radius S2- The three panels of the figure correspond to the three cases that need to be 
distinguished: (a) no overlap: r > + s 2 , (b) partial overlap: Si — s 2 < r < s 1 + s 2 , and (c) 
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complete overlap, r < s\ — s 2 . In all three panels of Fig. [2j the rectangle ABCD is the inte- 
gration domain, and the diagonal line through points E and H is the line of non-analyticities 
(where r — x — y = 0). For points below this line, the absolute value in the argument of 
nn in Eq. (27) may be omitted, while for points above this line, it changes the sign of the 
argument. Considering first case (a), i.e., no overlap, one sees from Fig. [2]^a) that 

Vij( r ) = -^AEH — -^BEG — -^DFH + -^CFG' (28) 



where /xyz * s the integral (27) with the absolute value sign omitted, and evaluated over the 
area of the triangle XYZ. For case (b), i.e., partial overlap, one finds from Fig. [2~ 



J AEH J BEG J DFH ' J < 



CFG' 



(29) 



where the superscript "— " indicates that the sign of the argument of 
changed. Finally for case (c), one finds from Fig. [2]^c) 



V i3 {r) 



AEH 



BEG 



DFH 



CFG- 



in Eq. (27) is 



(30) 



Note that for even basic potentials pn and nn , the sign of the arguments is inconsequential, 
so that all three cases (28)-(30) will have the same functional form. 

The integration limits appropriate for the triangular regions are easily determined from 
Fig. [2] This yields the following explicit expression for the integral /aeh : 



'AEH 



r+si 



r J-s 2 



dx / dy Ki(y, s x ) Kj(x, s 2 ) (r-x-y) nn (r -x-y) 



Aij(r, si,s 2 ). 



(31) 



Here, the identification with followed from Eqs. (20) and (22). Given the form of the 
auxiliary potential in Eq. (31), it is not hard to show that 



J BEG 



Aij{r,s^ -s 2 ) 



'DFH 



Aij(r, -si, s 2 ) 



'CFG 



Aij(r,-s h 



■82) 



(32) 



so that with Eq. (28) one finds for the non-overlapping case 



Vij(r) = Aij(r, s 1; s 2 ) - Aj{r, s 1? -s 2 ) - Aj{r, -si, s 2 ) + A {j (r, -s 1 , -s 2 ) 
= Ay(r, [si], [s 2 ]), 



(33) 



As was the case for Aij may have divergent parts which cancel in Eq. (33) and will be 
omitted below. 

According to Eqs. (28) and (29), the partially overlapping case (b) only requires replacing 

iven by 

dx / dyKi(y, si) Kj(x, s 2 ) (r-x-y) 4> n n(-r + x + y). (34) 



j cfg b y ^cfg> whicn is s iven b y 

1- - - r 

J CFG - / 

I Jl 51 

Substituting y — > —y, x - 
so that for d < r < D: 



x, and using that Ki and K~ are even in x and y, one finds 



-^CFG — Aij( 



-r,si,s 2 



(35) 



Vij(r) = Aj{r, si, s 2 ) - Aij(r, s u -s 2 ) - A„(r, -s 1 , s 2 ) + A^-(-r, s u s 2 ) 

= A ij ((r),s 1 ,s 2 ) - A ij (r,(s 1 ,-s 2 )) (36) 
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For the fully overlapping case, finally, one furthermore needs to replace /^ph by 



'dfh 



T Jr—si 

A 



si 



dx / dy Ki(y, si) Kj(x, s 2 ) (r - x - y) nn (-r + x + y) 



-r, si,-s 2 



whence for r < d: 



Vij(r) = Aij(r, si, s 2 ) - A^r, si, -s 2 ) - Ay(-r, s 1( -s 2 ) + Aij(-r, s 1( s 2 ) 



(37) 



(38) 



The reason that this is not symmetric in s\ and s 2 is because of the assumption that s± > s 2 . 
With si < s 2 and r < s 2 — si, one would have obtained V^-(r) = Ay((r), [si],s 2 ). This 
completes the derivation of Eq. (18). 



There is a degree of freedom in choosing the auxiliary potentials in Eqs. (17) and (18) 



since they enter only in specific combinations. In particular, according to Eq. (17), the 



effective point-nanoparticle potential is either r symmetric or s-antisymmetric. Thus, one 
may replace Ai(r, s) by Ai(r, s) + X(r, s) if the function X(r, s) is antisymmetric in r as well 
as symmetric in s, i.e., if 



X{r, s) = X{r, -s) = -X(-r, s). 



(39) 



Conversely, any terms in A4 that satisfy Eq. (|39[) are irrelevant to Eq. (17) and may, therefore, 



be omitted. Similarly, the effective inter-nanoparticle potential in Eq. (18) is not affected 



by adding a function Y(r, si, s 2 ) to the auxiliary potential A^, as long as Y satisfies 

Y(r, si, s 2 ) - Y(r, -s 1 , s 2 ) - Y(r, s 1 , -s 2 ) + Y(r, -s 1 , -s 2 ) = 

Y{r,s 1 ,S2) = Y{-r,-s u -s 2 ), (40) 

while terms present in A^ that satisfy these relations are irrelevant, and may be omitted. 



IV. SOLID AND HOLLOW NANOPARTICLES 

Two particular cases of the internal nanoparticle densities p will be considered in detail 
below. The first is a uniform internal density p inside a solid sphere of radius s: 



p(x) = pQ(s — x). 



(41) 



Since Eq. (41) is of the form 0,6(3 — x)x l with i = and a = p, Eq. ^ gives for the 
effective point-nanoparticle potential 



V pa (r) = pV (r). 



(42) 



Similarly, the effective inter-nanoparticle potential of two solid nanoparticles of uniform 
density p\ and p 2 , and radii s\ and s 2 , respectively, satisfies [cf. Eq. ([7|] 



V m (r) = Pip 2 V r 00 (r). 



(43) 



The second type of "internal" density p(x) considered here is that of hollow nanoparticles, 
whose density is concentrated on the surface of the sphere, i.e., 



p(x) =pd(s-x), 



(44) 
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where p is the surface density on the area of the sphere of size s. This density is appropriate 
to describe e.g. buckyballsP^ The density in Eq. (44) cannot be written in the form Eq. (|3]), 
but it is linked to the uniform internal density in Eq. (41) by 



pS(s — x) = p 



_ 8B(s — x) 



ds 



(45) 



Consequently, the effective point-nanoparticle potential for this case is given by 

V pn (r) =pV h (r) } (46) 

with 

= (47) 



ds 

where the subscript h indicates that this potential acts between a hollow nanoparticle and 
a point particle. 

In a similar fashion, the inter-nanoparticle potentials for a solid and a hollow nanoparticle 
(sh) is given by 

Kn(r) =Pip 2 V sh {r) (48) 
and the potential for two hollow nanoparticles (hh) satisfies 



V m (r) 



Pip 2 V hh (r), 



(49) 



where pi and p2 are the surface density of the two nanoparticles, while the scaled inter- 
nanoparticle potentials in Eqs. (48)-(49) are given by 

dV 00 (r) 



V sh (r) 
V hh (r) 



ds 2 
d 2 Voo(r) 

dsids 2 



(50) 



Thus, the effective potentials V^, V s h and Vhh can be found by differentiation once V and 
Voo, are known. 

The effective potentials for solid nanoparticles can be expressed in terms of auxiliary 
potentials A and A 00 using Eqs. (17) and (18). In applying Eqs. (47) and (50) to these 
expressions, it should be realized that taking a derivative turns an antisymmetrized function 
into a symmetrized one, and vice versa. Thus, by defining 



A h (r, s) 
A sh (r, s 1 ,s 2 ) 
A hh (r,s 1 ,s 2 ) 
one gets for the effective potentials 

T/ (r \ _ / 4 h ((r),s) if r < s 
Vh[r) ~ \ A h (r,(s)) i£r>a, 



9A (r, s) 
ds 

dA 00 (r, si,s 2 ) 
ds 2 

d 2 A 00 (r,si,s 2 ) 
ds 1 ds 2 



(51) 



(52) 
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V sh (r) 



V hh (r) 



Ah((r),[s 1 ],s 2 ) if r < \d\ and s 1 < s 2 

A S h((r), si, (s 2 )) if r < <i and si > s 2 

A sh ((r), si, s 2 ) + A sh (r, [s ls -s 2 ]) if |d| < r < L> 

A sh (r, [si],(s 2 )) ifr>D, 

Ahh{{r), (si), s 2 ) if r < |d| and s a < s 2 

Ahh((r),si,(s 2 )) if r < and si > s 2 

+ A hh {r,{ Sl ,-s 2 )) if |d| < r < £> 

Ah(r,(si),(s 2 )) ifr>D. 



(53) 



(54) 



V. EFFECTIVE POTENTIALS FOR UNIFORMLY SOLID AND HOLLOW 

NANOPARTICLES 



Power laws 



Pair potentials of power law form 



(55) 



with n integer, are basic building blocks of many atomic and molecular pair potentials, such 
as the Coulomb potential (n = 1) and the Lennard- Jones potential (a linear combination of 
n = 6 and n = 12). Note that here and below, a superscript on a potential represents an 
index, not a power. 

The effective potential Vq for a point particle and a solid nanoparticle of radius s whose 
points interact with the particle through pn = <p n is given in terms of the auxiliary potential 



by Eq. (17). The auxiliary potential follows from Eqs. (19), giving, for general n, 

s 2 — (r — y) 2 



A n (r,s) 



71 



r Jo 



r+s 



dy 



y 



n-l 



2vr[r + (n- 3)s] 
(n - 2)(n - 3)(n - 4) r (r + s) n - 3 ' 



where divergent terms were omitted, as explained in Sec. Ill 



(56) 



2, 3 



The right hand side of Eq. (56) becomes ill-defined for the specific values n 
and 4. This is caused by a term proportional to x 11 ~ n ~ 1 in the integrand in Eq. (56) (with 



n' = 2, 3 or 4), which when n = n' should have resulted in a term ln(r + s) instead 
of the erroneous and ill-defined expression ^±^- 

[{n — n') x "_ " ] = lnx, this can be fixed by substituting 



that occurs in Eq. (56). Using that 



lim„_ >n . ^ 



A n ' — ► lim ^-\(n-n')A r 
n^n> dn 



(57) 



Applied to Eq. (56), this gives 



A 2 (r,s) 



A 3 {r,s) 



At(r,s) 



rx{r + s)(3r — s) 



2r 



+ 



tt(s 2 — r 2 ) 



ln(r + s) 



2tts 



+ 2vr ln(r + s) 



r 

7r(3r + s) it 

■—. f ln(r + s) 

2r(r + s) r 



(58) 
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The effective potential V$ is obtained from these expressions for the auxiliary potential 
using Eq. (17). 

From Eqs. (51) and (56), it follows that the auxiliary potential for a hollow nanoparticle 

2ns 



and a point particle is given by 



Al{r,s) 



(n-2)r(r + s) n ~ 2 ' 



Equation (59) is ill-defined for n = 2, in which case one uses Eq. (57) to find 



2>KS 



ln(r + s). 



(59) 



(60) 



The effective potential V£ is now obtained from Eq. (52). 



For the effective inter-nanoparticle potential Vqo, the auxiliary potential formulation (18) 



holds with i — j — 0, where the auxiliary potential is found using Eq. (20) with 
giving 



47T 2 p n (r, si, Sg) 

(n — 7)(n — 6)(n — 5)(n — 4){n — 3) {n — 2) r (r + s\ + s 2 ] 



n-5 ' 



where 



p n (r, s 1 ,s 2 ) = r 2 + (n - 5)(s 1 + s 2 )r + (n - 6)[sl + si + (n - 5)sis 2 ]- 



(61) 



(62) 



The expression in Eq. (61) is ill-defined for n = 2, 3, 4, 5, 6 and 7. Using again Eq. (57), 
the correct expression for A^ for these values of n is found to be 



4vr 2 



r (r + si + s 2 ) n ~ 5 H\ = 2 (l-n) 

r 7 1 

Pn(r, si, s 2 ) ln(r + s x + s 2 ) - JZ~ 

1 = 2 1 n 



~ s i - s i - ( s i + s 2) r + (11 _ 2n)s 1 s 2 



(63) 



According to Eq. (51), the auxiliary potential for a solid sphere of radius si and a hollow 
sphere of radius s 2 can be found by taking the derivative with respect to s 2 , yielding, for 
general n, 



-An 2 s 2 [r + {n- A)s 1 + s 2 ] 



(64) 



(n - 5)(n- 4)(n - 3)(n - 2)r(r + si + s 2 ) n_4 ' 

Finally, the effective potential for two hollow spheres follows from another derivative with 
respect to si [cf. Eq. (51)], leading to 

A£ h (r, si,s 2 ) 



47T 2 SiS 2 



(n - 3)(n - 2) r (r + si + s 2 ) n - 3 ' 



(65) 



For the ill-defined cases of Eqs. (64) and (65), one can use Eq. (57) to get expressions similar 
to the one in Eq. (63). 
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FIG. 3: Typical example of effective potentials based on an exponential interaction [Eqs. (18) and 
(69)]. The left panel shows the point-nanoparticle potentials for solid (s) and hollow (h) spheres 
with radius s = 3, while the right panel shows the inter-nanoparticle potentials for radii s± = 4 
and s 2 = 1. 



B. Exponentials 

The effective interactions as a result of the exponential pair potential 

E (r) = e~ r 



(66) 



will now be derived. Substituting this potential for pn in the expression ( 19 ) for the auxiliary 
potential gives 



, E . , 2vr(3 + r + sr + s 2 + 3s) _ r _ s A 
A*(r, s) = —± >-e r s + 47r, 



(67) 



where an irrelevant expression satisfying Eq. (39) was omitted. From Eqs. (51) and (67), 
the auxiliary potential for a point particle and a hollow nanoparticle is found to be 

2ns(l + r + s) „ „ 



(68) 



Note that the corresponding effective potentials follow from Eqs. (17) and (52). 



The effective inter-nanoparticle potential is of the auxiliary potential form (18) with 

Vn = E , giving 



j = 0. The auxiliary potential A^ Q is found using Eq. (20) with 



A E (r,s 1; s 2 ) =4tt 2 



(r + sj + s 2 + 5)(si + l)(s 2 + !) + !- s lS2 _ r _ Si _ 



+ 



Ti- 

3r 



r 

,2 i 2 



si + s 2 )(si + s 2 2 - Sl s 2 )r + 6(sf + - 4)(r 2 + 4) - r 4 + 3(sJ - s 2 2 ) 2 + 24 



(69) 



where an expression satisfying Eq. (40) has been omitted. Using Eqs. (51), the auxiliary 
potential for the interaction between a solid and a hollow nanoparticle and between two 
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FIG. 4: Typical example of effective potentials based on the London-van der Waals interaction, 



i.e., the power law in Eq. (55) with n = 6. The left panel shows the potential for a point particle 



and solid or hollow nanoparticle of radius s = 3, the right panel shows the potential for two 
nanoparticles of radius s\ = 4 and s 2 = 1. 



hollow particles are found to be 



-4n 2 s 2 [(r + si + s 2 + 4)(si + 1) - Sl ] 



r 

| 47r 2 s 2 [(r + s 2 ) 2 -s 2 + 4] 
r 

47r 2 sis 2 (r + si + s 2 + 2) 



8tt 2 s 1 s 2 



(70) 
(71) 



Figure |3] shows a typical example of the effective potentials derived from the exponential 



basic potential [cf. Eqs. (17), (18), (52)-(54) and (67)-(71)]. One sees that these effective po- 
tentials are very smooth and do not have a hard core, which is typical for effective potentials 
based on a basic pair potential that does not diverge for small distances. 



C. Examples using common pair potentials 

London-van der Waals potential 
In this section, the effective potentials based on the London-van der Waals potential 



<f{r) 



(72) 



will be presented. Note that the negative prefactor that occurs in front of the attractive 



London- van der Waals interaction has been omitted here. Substituting n = 6 into Eq. (56) 



and using Eq. (17), one finds the London- van der Waals potential for a solid nanoparticle 



and a point particle: 



47TS 3 
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for r > s. This effective potential becomes infinite for r < s. For the London- van der Waals 
interaction of a hollow nanoparticle with a point particle, Eqs. (52) and (59) with n = 6, 
lead to 

Airs 2 8vrs 4 



(r 2 — s 2 ) 3 



(r 2 -s 2 Y' 



(74) 



The effective London-van der Waals interaction potential for two solid nanoparticles is 
determined by substituting n = 6 into Eq. (63), and using Eq. (18), which gives 

(75) 



7T SiS 2 



7T SiS 2 



3(r 2 -d 2 ) 3(r 2 
This result coincides with that of HamakerP^ 



D 2 



7r r — JJ 
6 r 2 — d 2 



Using Eqs. (50) and (75), or using Eqs. (64) and (53), one finds for the London- van 
der Waals potential for a solid nanoparticle of radius s\ and a hollow nanoparticle of 
radius s 2 

2tt 2 sis 2 D 2n 2 s 1 s 2 d 



7T 2 S 2 



7T 2 S 2 



The effective London-van der Waals potential for two hollow nanoparticles, finally, 
is obtained from Eq. (76) using Eq. (50), or alternatively from Eqs. (65) and (54), with the 
result 

6 87r 2 gig 2 D 2 _ 87r 2 S!S 2 d 2 2ir 2 s 1 s 2 _ 2tt 2 SiS 2 

hhlH " - 3 ( r 2 _ £ 2 ) 3 3 ( r 2 _ rf 2)3 + 3 ( r 2 _ £,2)2 ^2 _ rf 2)2 ' \ > 

Figure [4] shows a typical example of the effective potentials for the London-van der Waals 
interaction as the basic pair potential. 



Morse potential 

The Morse potential 



1M1 



r ) =e -26(r-l)_ 2e -6(r-l) ) (7g) 

is used e.g. for molecular bonds and for pure metals.^ It is a sum of two exponential 
functions, so having derived the formulas for the exponential potential in Sec. VB[ one 



easily finds the corresponding point-nanoparticle interactions by taking the combinations 

C(0 = V*(2br, 2bs) - -A V E (br, bs) (79) 

2b o 6 

V^(r) = ^ K E (2&r, 2bs) - — V h E (6r, te), (80) 
where the notation V E (ar, (3s) indicates that in V E and V E , r is to be replaced by ar and s 



by (3s. Likewise, the inter-nanoparticle interactions for the Morse potential in Eq. (78) are 
given by 

e 2b 2e h 
V «>( r ) = ^ K>o(2&r, 2b Sl , 2bs 2 ) - — V E (br, b Sl , bs 2 ). (81) 

2b 2e b 
(0 = ^^(2^,2^,26^) - — V£{br, b Sl , bs 2 ). (82) 

2b p b 

V&{r) = ^ V*(2br, 2bs u 2bs 2 ) - — V*{br, bs 1 , bs 2 ). (83) 
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FIG. 5: Example of Morse effective potentials for b = 2.6. The left panel shows the effective 
potential for a particle and a solid or hollow nanoparticle of radius s = 3, the right panel shows 
the effective potentials for two nanoparticles of radius s± = 4 and S2 = 1. 




1 2 3 4 5 6012345678 



r r 

FIG. 6: Example of Morse effective potentials for b = 5.6, for which the Morse potential resembles 
the Lennard-Jones potential. The left panel shows the effective point-nanoparticle potentials for 
s = 3, the right panel shows the effective potentials for two nanoparticles of radius s\ = 4 and 
S2 = 1. Note that V s h and Vh s are nearly the same for r > D. 



Two examples of the Morse-based effective potentials are shown in Figs. [5] and |6j for 
b = 2.6 and b = 5.6, respectively. For the lower value of b, there is a low barrier for a point 
particle to penetrate a nanoparticle as well as for one nanoparticle to penetrate another 
(cf. Fig. [5|, while for the larger value of b this is virtually impossible (cf. Fig. pj if the 
energies of the particles are of order 1. 
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Buckingham potential 
The modified Buckingham potential 



0B (r) = J° o 6 [t f r<r : (84) 

r v 1 ae br — cr b if r > r* v 7 



is made up of an exponential part, for which the results of Sec. VB apply, and an attractive 
London-van der Waals term treated above. In addition, one needs to take the cut-off r* into 
account. This cut-off is necessary because otherwise, for small enough r, the Buckingham 
potential would become negative. Thus, the effective point-nanoparticle potentials are 

\/B/ \ _ / oo if r < s + r* , , 

[T) ~ \ §V*{br, bs) - cV 6 {r) if r > s + r* ^ } 

[ §2V^(br, bs) - cV£(r) if r < s - r* 
V*(r) = I oo if \s-r\ < r* (86) 

1 £ Vh( br , bs ) - cV h( r ) if r > s + r* 

while the effective inter-nanoparticle potentials are given by 

v B (r) _ / if r <D + r* 

001 J \ f V^(6r,6si,6s 2 ) - C y 6 (r) otherwise 1 J 

t/B.x _ / oo if -d - r* < r < D + r* 

shl ' \ £Vi(br,b Sl ,bs 2 ) - cVl(r) otherwise { } 

r/B/x/oo if |d| -r* <r<D + r* 

hhl J \ ^V^(6r,6 Sl ,6s 2 ) - C y h 6 h (r) otherwise. 1 ' 

While the effective potentials due to the exponential pair potential are different for different 
cases (no overlap, partial overlap, and complete overlap), because of the presence of a cut-off 
r*, only the non-overlapping case is relevant here. 

In Fig. [7], a typical example of these potentials is shown. Note that while it is possible 
for a point or nanoparticle particle to be inside the hollow nanoparticle (as long as there is 
no overlap), there is an infinite barrier to get inside from the outside, in contrast with the 
effective potentials based on the Morse potential. 



Lennard-Jones potential 



One of the most often used potentials in molecular dynamics simulations is the Lennard- 
Jones potential,^ which in reduced units reads 



1 

^12 



(j) 12 (r) - 2(j) b (r). 



(90) 



Since the attractive part of the Lennard-Jones potential in Eq. (90) was handled above, one 

12 to find the effective potentials for Lennard-Jones 



only needs to add the repulsive part r 
nanoparticles. Substituting n — 12 into the results of Sec. |V A 



and using the relations 
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FIG. 7: Typical example of effective potentials based on the Buckingham potential for a = e' 
b = 13 and c = 2, with the cut-off r* set to 1/4. The left panel shows the effective potential 
for a point particle and a solid or hollow nanoparticle of radius s = 3, the right panel shows the 
potentials for two nanoparticles of radius s\ = 4 and s% = 1. 



between auxiliary and effective potentials, one finds 



K 12 (r) 



47TS 3 



3(r 2 - s 2 ) 

47TS 2 



+ 



80vrs 9 + 432vrrV 



(r 2 - s 2 ) 6 



+ 



45(r 2 - s 2 ) 9 

6, 
5' 



64vrrV(r 4 + |s 2 r 2 + s 4 



(r 2 - s 2 ) 10 



(91) 
(92) 



7T 



37800r 



(r + ?£) 2 + f £ 2 - f ci 2 
(r + L>) 7 

5 n 2 15 j2 



D) 



2 + |£> 2 - f d 2 



(r + \df + f rf 2 - f D 2 

(r + rf) 7 
(r - |rf) 2 + f d 2 - fD 2 



TT 2 S 2 

1260r 

2tt 2 SiS 2 
45r 



(r-D) 7 
r+\D + ld r-\D 



7 



(r - d) 7 

9j i 7 



(r + £>) f 



(r - D) 8 (r + d) 8 

1111 
(r + D) 9 + (r-D) 9 ~ (r + d) 9 ~ (r - d) 9 



^ + r + ^d + \D r-ld-\D 



(r - dy 



(93) 
(94) 
(95) 



The potential Vqq is in agreement with the result in the appendix of Ref. [U 

The point-nanoparticle potentials for the Lennard- Jones potential are now given by 



V 12 (r) - 2V 6 (r) 



4ns 3 



3(r 2 - s 2 ) 



+ 



80ns 9 + 4327rrV 



57TS" 



45(r 2 - s 2 ) 9 3(r 2 - s 2 ) 3 



V£ 3 (r) = V^(r) - 2V°(r) 



Ans z 



(r 2 — s 2 ) 



+ 



64rrr 2 s 4 (r 4 + |sV + s 
(r 2 - s 2 ) 10 



ins 



lQns 



-^2 q2 j3 



fr 2 -s 2 ) 4 ' 



(96) 



(97) 



Equation (96) is a more concise notation of the result of Roth and Balasubramanya [Eq. (2) 
in Ref. H3] ■ Likewise, the inter-nanoparticle interactions due to a Lennard- Jones potential 
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FIG. 8: Typical effective potentials based on the Lennard- Jones potential [Eqs. (96)-(98)]. One 
the left, the potential for a point particle and solid or hollow nanoparticle of radius s = 3 is shown, 
and on the right, the potentials for two nanoparticles of radius s% = 4 and S2 = 1. 



U_ 1 2 3 4 5 678 

ss 2" 219136 24064 3456 _2^ 2f 

"i 315 4725 675 225 45 9 

«sh 2^ _2^ 2^£ _2_^ 2j_ _2^ 2i 

"« 315 315 45 45 3 15 3 

„,hh 2^ _2^ 2^ 917504 57344 13312 14336 9 7 o4 

"i 45 5 5 30 5 5 15 z z 

TABLE I: Coefficients for the polynomials appearing in the effective inter-nanoparticle potentials 
based on the Lennard-Jones potentials, i.e., Vqq 1 , Ksh an< ^ ^hh J m ^qs. (99)-(101). 



are given by 

V^{r)=V l f{r)-2V^r\ (98) 

where ij = 00, sh or hh. In Fig. |8j a typical example of these effective potentials is shown. 
Note the hard core part of the potentials. For the specific case of system of nanoparticles with 
the same radii si = S2 = s, studied in Ref. [23J, the effective inter-nanoparticle interactions 
can be written in terms of rj = r/s as 



ss „2i 



s e v s (^2 _ 4 )7 
S 7 T] 8 {jf — 4) 8 



4vr z 



V 2 -2 



7T 



7T 2 E 8 



i=0 



af 1 r] 2i 



3 T] 2 (r] 2 - 4) 

32tt 2 
3s r/ 2 (r? 2 — 4) 2 
2 r^ 4 + 6?? 2 



In 1 



4 



(99) 
(100) 

4)9 — S 2 V 4^2 _ 4 )3 ( 101 ) 

with the a coefficients given in Table [ij Equation (101) is the so-called Girifalco potential.^ 



s 8 rj 10 (r] 2 



- 32rr 



VI. ACCURACY OF THE LENNARD-JONES BASED EFFECTIVE 
POTENTIALS FOR FCC NANOPARTICLES 

Since the effective potentials derived above are intended to model nanoparticles, it is 
natural to ask to what extent they can represent the interactions of nanoclusters composed 
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FIG. 9: Comparison of the fitted radius s and the a priori radius s* of the fee nanoparticles. The 
fit is based on minimizing A pn (s), but minimizing A nn (s) instead gives indistinguishable results. 



of atoms. This obviously will depend on the structure of the nanoclusters, but to get at 
least a partial answer, the fec-based nanoparticles of Sec. [IT] will be used again, with the 
basic pair potentials <f) pn and <ft nn given by LJ in Eq. (90). This potential has a minimum 
at r = 1, which sets the unit of length. The fee nanoparticles are constructed from an fee 
lattice with mean density p = 1 by picking an atom and including all atoms within a given 
distance from it. Note that this gives only specific values for the number M of included 
atoms, since many atoms lie at the same distance in the crystal structure. Here, M will be 
restricted to less than 20,000, resulting in 206 clusters, the largest of which has M =19,861 
atoms. 

The mean density p = 1 for the fee nanoparticles is not unrealistic: It results in a lattice 
distance a = 4 1 / 3 (Ref. UHl p. 12), i.e., the ratio of the lattice distance to the interaction 
range is 4 1 / 3 rs 1.587. This is comparable to the case of platinum nanoparticles in water: 
Assuming the lattice distance a is the same as in a bulk platinum crystal, a = 3.92 A 
(Ref. UHl p. 23), and using that the interaction range of Pt atoms with water is of the order 
of 2 to 3 Aj 23 one finds a similar ratio of 3.92A/2.5A = 1.568. 

To test the applicability of describing these fee nanoclusters as spheres with a constant 
density, one should compare the effective point-nanoparticle potential V pn = pV^ 3 to the 
result of summing the potentials </> LJ between the point particle and each of the atoms in the 
fee nanoparticle. Similarly, the effective potential V nn = p 2 V^ between two equally sized 
nanoparticles should be compared to the result of summing the potentials between the each 
of the atoms of one of the nanoparticles with each of the atoms in the other. 
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FIG. 10: Deviations A pn and A nn of the effective potentials from the atom-by-atom summed 
potentials for the fee nanoparticles as a function of the fitted radius s. 



However, there are two difficulties in performing these comparisons. First, the effective 
potentials are spherically symmetric, but the summed potentials will not be, since the fee 
nanoparticles are not truly spherically symmetric. Therefore, the comparison will be made 
with the summed potentials averaged over all orientations of the nanoparticles, which will 
be denoted by V™ m and V^ m . 

The second problem with the comparison is that the radius s of the nanoparticle, which 
is a parameter in the effective potentials, is not well defined. A reasonable a priori radius 
would be s* = [3M/(47rp)] 1//3 , but other values for the radius s close to s* are just as 
reasonable. Thus, the radius may be viewed as a fitting parameter, which will be adjusted 
to minimize the difference between the effective and the summed potential. To be precise, 
the following quantities are minimized by varying s: 



A pn = II dr 



A nn = < I dr 



V:r(r)- P 2 (s)V y(r)-\ (102) 



1/2 



1/2 



Here, p(s) = 3M/(47rs 3 ), and the prime denotes the restriction on the integration that 
Vp™ m (r) < 3V* Q or V^ m (r) < 3V* n , respectively, where V* Q and V* n are the absolute value 
of the minima of V^ m and V^ m . The restriction is needed to make the integrals converge. 
The results depend very little on the precise choice of the restriction. For instance, changing 
the restriction to 2V* instead of 3V*, shifts the values for the radii s only by an amount of 
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FIG. 11: An example of very good agreement between the effective and summed potentials, which 
occurs for an fee nanocluster of size M = 18053 with an effective radius of 16.27 (in dimensionless 
units). Crosses represent the orientationally averaged summed potentials V p ^ m (left) and V r 

im /'»() ,lllu •'nn = P 2 ^00 



sum 
nn 



(right), while the solid lines are the effective potentials V pn = pV^ (left) and V nn = p 2 Vqq (right). 



the order of 1CT 4 . 

The values of the radius that result from minimizing A pn for the 206 cluster configurations 
with M <20,000 are shown in Fig. [9] It is seen that with the exception of some of the 
smaller clusters, the values of fitted radii s typically lie close very to the a priori radius s*. 
Minimizing A nn instead results in the same values for the radii to within 0.3%. 

To get an idea of the accuracy of the fit as a function of the size of the nanoparticles, one 
may investigate the values of the dimensionless deviations 

a _ ^pn , a _ A nn 



Pn R 1/2 V* ' E 1/2 V* ' 

-fi-pn Vp n -Ttnn V nn 

The length scales R pn and R nn are chosen as the lengths of the intervals contributing 99.9% of 
the values of the integrals in Eqs. (102). This typically gives R pn ~ 1.35 and R nn ~ 2 for the 
size of clusters investigated here, and these values of R pn and i? nn were used for all clusters. 



The dimensionless deviations are plotted in Fig. 10 One sees a high degree of correlation 
between the accuracy of the effective potential for a nanoparticle and point particle and 
the accuracy of the effective inter-nanoparticle potential. The deviations are furthermore 
typically small, indicating that their is good agreement between the effective potentials and 
the sum of atom-atom potentials, although the deviations are larger for specific cluster sizes. 



As extreme examples, Fig. 11 shows a case of very good agreement and Fig. 12 shows a case 
of poor agreement. In these figures, the effective potentials and the summed potentials are 
compared for M = 18053 with s = 16.27 and M = 17357 with s = 16.04, respectively. 
Note that the agreement is never very bad, but for the latter, the depth of the minimum is 



somewhat underestimated by the effective potentials, as the insets of Fig. 12 show. 

It is hard to say in general why the smooth, constant density description works better for 
some clusters than for others. For some of the smaller nanoclusters with poorer agreement, 
inspecting the spatial structure of the nanocluster shows a rather rough surface, which could 
be the explanation. But for the larger nanoparticles, such a difference in roughness is hard 
to distinguish. 
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FIG. 12: An example of poorer agreement between the effective and summed potentials, which 
occurs for an fee nanocluster of size M = 17357 with an effective radius of 16.04 (dimensionless 
units). Solid lines represent the effective potentials V^n = pV^ (left) and V nQ = p 2 V g J (right), 
while crosses are the orientationally averaged summed potentials y p s ™ (left) and V^ m (right) 
which result from the sum over the atoms. The inset in the right plot zooms in on the minimum, 
and shows that its depth is underestimated by the effective potential. 



VII. DISCUSSION 



A general effective description for nanoparticles was presented, starting from a smoothing 
procedure in which the real spatial density profile inside the nanoparticles is replaced by a 
spherically symmetric one. The resulting effective interactions between a nanoparticle and a 
point particle as well as between two nanoparticles are then given by spherically symmetric 
potentials, thus greatly simplifying the description over an all-atom model. 

The main results of this approach are the formulation of the effective potentials in terms 
of auxiliary potentials, Eqs. (17) and (18), which provide a unified description of overlap- 
ping and non-overlapping configurations. The auxiliary potentials are related to the basic 



interaction potentials through Eqs. (19) and (20). Furthermore, the effective potentials for 



hollow particles were found to be related to those for solid nanoparticles by simple differ- 



entiation with respect to the radii of the nanoparticles, see Eqs. (47) and (50), and as such 



also allow a formulation in terms of auxiliary potentials, as given in Sec. IV 



As an application of the formalism, explicit effective pair potentials for solid and hol- 
low nanoparticles were obtained for various basic pair potentials. Different pair potentials 
have different applications. For instance, the Lennard- Jones potential is a general-purpose 
potential, while the Buckingham potential is suited to describe the physics of particles 
close together such as in high pressure systems. These basic potentials result in effective 
nanoparticle potentials with hard cores plus a soft potential. They reduce in limiting cases 
to some of the existing model potentials for colloids, such as hard spheres and the Hamaker 
potential,™™ 1 but not to more ad hoc models such as the description of a colloid as a single 
big Lennard- Jones particle.' 22 ' In contrast, the Morse potential is able to describe bounded 
systems or penetrable particles, making it possible to model nanoparticles that could pas- 
sively capture and trap specific types of particles. This could have applications in modeling 
drug delivery by nanoparticles^ and viral capsidsP 

For the case of a Lennard- Jones basic potential, a comparison was carried out with 



25 



an atomic model of a nanocluster. In this model, the atoms making up the nanoparticle 
were assumed to be arranged in an fee lattice structure. To find approximate spherical 
structures, the atoms were restricted to lie within a certain distance from the central atom 
in the nanocluster. Configurations with up to 19,861 atoms were studied. The effective 
potentials were compared with the orientionally averaged sum of Lennard- Jones potentials 
due to the individual atoms. The agreement tends to be very good, provided the radius in 
the effective description is treated as a fitting parameter. For some configurations, however, 
the fitting procedure underestimates the depth of the minimum of the potentials. This may 
be due to surface roughness of these structures, which is caused by the imposed fee structure 
and unlikely to be relevant for real nanoclusters. 

The application of the explicit expressions for the effective potentials to numerical studies 
of spherical nanoparticles is in principle straightforward. In fact, the potentials in Eqs. (96) 
and (99) have already been used in a numerical study of single particle transport in an 
equilibrium nanofluid composed of solid nanoparticles and fluid particles interaction through 
Lennard- Jones interactions, where the validity of a Gaussian approximation of the Van Hove 
self-correlation function was investigated, and found to hold up to picosecond time scales 
for the fluid particles, and up to five to ten times longer (depending on temperature) for 
nanoparticles with a size of about 2 nmP 

Given the explicit expressions for the effective potentials, the description allows a fairly 
direct route toward a qualitative model for a given system of nanoparticles in a fluid, since 
reasonable values for the parameters for commonly used pair potentials are available in the 
literature,^ while the number of atoms in a nanoparticle and its radius could be taken 
from experiments or theoretical calculations.^ Furthermore, the effective potentials have a 
physical range based on the interaction of their constituents rather than on their radius. 
Therefore, the effective potentials that were derived here are expected to be useful for the 
qualitative description of a wide variety of systems, from mono-disperse nanoparticles in a 
fluid to mixtures of different kinds of fluid particles, nanoclusters or buckyballs. 

A number of interesting extensions present themselves for future research. For instance, 
while the nanoparticles were assumed to be composed of one kind of particle only, poten- 
tials for nanoparticles composed of several types of particles can also be derived within the 
current context if the distribution of the types is either homogeneously mixed or distributed 
in spherical shells (so-called core-shell nanoparticles^^). The spherical symmetry of the 
effective potentials, which decouples the rotational and translational degrees of freedom, 
could be lifted to extend the model to include rotational motion. This may be done by 
adding interaction sites on the surface of the nanoparticle or a multipole expansion. As long 
as the orientationally dependent potential is available, there are no obstacles in molecular 
dynamics simulations of such systems.^ Furthermore, combining the current model with the 
mesoscopic fluid model of Malevanets and KapraP^ would yield a numerically efficient model 
of larger nanoparticles and colloids that includes hydrodynamic effects. These avenues are 
currently being investigated. 
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APPENDIX A: THE KERNEL K, 



■i j 



The integral in the expression for the kernel Ky in Eq. (14) will be worked out now. 
Using Eq. ^ and the binomial formula for (x — y) l+2 , one finds, after resummation, that 
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(A2) 



where yt = max(— s 2 ,x — st) and y 2 = min(s2,a; + st), which are due to the finite support 
of the kernels Ki and Kj, and F is the hypergeometric function.^ Despite its complicated 



appearance, Eq. (A2) is simply a piecewise polynomial in x of degree i + j ' + 5 at most. To 
see this, it is useful to distinguish the following four non-trivial cases: case 1: x > and 
\d\ < \x\ < D, for which y x = x — S\ and y 2 = s 2 ; case 2: d > and |x| < \d\, giving yt = —s 2 
and y 2 = s 2 ; case 3: d < and \x\ < \d\, giving yt = x — St and y 2 = x + St] and case 4 : 
x < and \d\ < \x\ < D, for which yt = —s 2 and y 2 = x + st- There are in fact only two 
independent cases, because case 3 can be obtained from the result of case 2 by interchanging 
st and s 2 as well as i and j (which will also flip the sign of d), while the result for case 4 
can be obtained from that of case 1 by setting s x to — s 2 , s 2 to — st and introducing a minus 
sign, as can be proved by changing the integration variable in Eq. (Al) from y to x — y. 



Thus, one only needs to consider the cases 1 and 2. Changing the integration variable from 
y to z = s 2 — y and using the binomial formula, Eq. (Al) for case 1 yields 



j—m+k 

4 



( f) _ T \m+3 min(i+2,m+l) /_ \j+2-fe 

K ij (x,s 1 ,s 2 ) = (i + 2)\(j + 2)\ £ [U £ " 2 ^,„ (A3) 



which is polynomial in x of degree i + j + 5, while for case 2 the integral in Eq. (Al ) can be 
found by using the binomial formula for (x — y) l+2 , giving a polynomial of degree i + 2, i.e. 



K ij (x,s 1 ,s 2 ) 
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